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1 . SCOPE 


This specification establishes the modifications to the CLASY program as 
specified in lOSD Category 1 Job Order 81-127 Task Agreement, titled CLASY 
Program Modification. 



/ 


2. APPLICABLE DOCUMENTS 


The following documents form a part of this specification: 

J. 0. 81-127 Task Agreement titled; CLASY Program Modification 
TIRE 76-0076 

Memorandum dated March 17, 1976, from Rice University, Institute of Com 
puter Sciences (D. L. Van Rooy) to Ken Baker/TF3, NASA- JSC; Reference: 
documentation of SUPER-SCRAM 


3. SYSTEM DESCRIPTION 


3.1 HARDWARE DESCRIPTION 


Not applicable— see uection 4 for the operational configuration required by 
the modified CLASY prog?*o!^ 

3.2 SOFTWARE DESCRIPTION 


The CLASY system of subprograms was originated by Dr. Michael Rassbach, a 
post-doctoral contractor for NASA-JSC, Earth Observations Division (TF). 

The purpose for his development of CLASY was the implementation of an itera- 
tive statistical clustering algorithm which had theoretical promise for 
application to classification of earth resources (image) data acquired from 
the LANDSAT (formerly ERTS) satellite. 

These modifications to CLASY are (1) to improve the program execution 
time, especially input/output overhead, in order to make it feasible to 
test and evaluate the CLASY program, and (2) to implement a randomized 
input data scrambling technique which was obtained from Rice University by 
the Earth Observations Division's Research, Test and Evaluation (RT&E) 

Branch (TF3). The scrambling of input data vectors is a necessity in CLASY 
due to the algorithm's sensitivity to correlated data values. 

The CLASY system of subprograms consists of the main driver program, CLASY, 
and 56 subprograms, not including the Univac systerri routines utilized by 
the program. Four of the subprograms (LOCK, GET, FREE, and BYTRAN) were 
originally programmed in Univac assembly language (SLEUTH II), the remaining 
subprograms and the main program, CLASY, are converted to the Univac Fortran V 
language from the former (original) Univac reentrant Fortran language, RFOR. 



The driver program for the clustering system is CLASY. The data handling 
subprogram for the system is CLASYl (with an alternate entry, CLASY2). 

CLASYl is reprogrammed to implement the required data input improvement and 
to implement the randomized data scrambling technique. CLASYl provides the 
data setup for acquisition on demand by the iterative statistical subpro- 
gram, STATIS. STATIS obtains the data from CLASYl via the alternate entry 
into the CLASYl subprogram, which is CLASY2. STATIS initiates the clustering 
procedure, operating on one pixel (data vector) at a time in setting up 
clusters and making the cluster split/combine decisions. Each pixel is 
examined 10 times by STATIS during the clustering proceoure. 

To implement the required modifications and the modifications suggested by 
the program originator to improve the program's reliability, the following 
routines in the CLASY clustering system received changes: CLASY, CLASYl, 
ADJUST, LOCK (formerly LOC), CLPR, CLDUMP, MISH, SETUP9, CMBKIO, ELIM, 

MULTI, SEPER, TR, STATIS, CLUST. 

The modified subprograms are discussed in the order of their use by the 
CLASY system. 

The overall CLASY system is flowcharted in Appendix A. Listings of the 
modified routines are shown in Appendix B. Sample output from the CLASY 
system is shown in Appendix C. 

3.2.1 SOFTWARE COMPONENT NG. 1 (CLASY) 

3. 2. 1.1 Linkages 


CLASY is the driver program of the CLASY clustering system. CLASY calls 
SETUP9, READTP, CLASYl, MULTI, and CLUMSMP. 


3. 2. 1.2 Inteifaces 


The common blocks INFORM, CLUSTR, CLUS, MISC, and STPAR and calling argu- 
ments are used 'in the program CLASY as interfaces with other routines in 
the clustering system. 

3. 2. 1.3 Inputs 

CLASY calls SETUP9, which reads the input supervisor (control) cards. The 
supervisor cards and their functions are described in the discussion of 
SETUP9 (section 3.2.2), and on the description of CLASY input (section 4.2.1). 

3. 2. 1.4 Outputs 


None. 


3. 2. 1.5 Storage Requirement 

Storage used: Code = 146g Data = 47047g 

3. 2. 1.6 Description 

CLASY is the driver program for the clustering routines. It was rewritten to 
(1) enable CLASYl to retrieve large blocks of data (dimensioned ARRAY (20 000)) 
from drum (2) to pre-calculate the amount of data RREAD (drum read sub- 
routine) will read from drum, and (3) to make DATAB (the array containing 
the scrambled data) and PV (the array passed to the clustering routine, 

STATIS) to reside in the same locations as the large data array, "ARRAY". 

All other logic in the program remains the same as in the pre-modified 
version of CLASY. 

3. 2. 1.7 Flowchart 



See Appendix A. 


3. 2. 1.8 


See Appendix B for program. 

3.2.2 SOFTWARE COMPONENT NO. 2 (SETUP9) 


3. 2. 2.1 Linkage 


SETUP9 is called from CLASY. SETUP9 calls NXTCKR and NUMBER, which are 
entry points in subroutine FIND. 


3. 2. 2. 2 


Interface 


j 

'j 

j 


Interface is accomplished through calling arguments and the following 
common blocks: INFORM, SUPCUM, and CLUSTR. 


3. 2. 2. 3 Inputs 
See Appendix C. 

3. 2. 2. 4 Output 

SETUP9 prints out a summary of the input to CLASY, and also prints an 
error message on the line printer if an invalid input card is detected. If 
an error is detected, SETUPS prints the following message "INVALID INPUT 
CARD— IGNORED" , processing continues. 


3. 2. 2. 5 Storage Requirement ^ 

i 


Data = 157g 


' 



Storage used: Code = 275p 


3. 2. 2. 6 Description 


SETUP9 reads and analyzes all cards input to the CLASY program. SETUP 9 was 
modified to add "0" (zero) as a new symbol to the symbol array “SMBLS". 

The following control cards are input to the modified CLASY program, to be 
analysed by SETUP9. In all cards, the "keyworn' begins in card column 1, 
and any parameters on the card are placed from card columns 11 through 72, 
inclusive, 

1. "CHANNEL" CARD (i.e., "CHANNEL 1,5,9,13") 

The "CHANNEL" card specifies the channel numbers to be used in clustering the 
mul ti -channel data vectors. At present the maximum number of channels 
allowed to specify is fr-^r. The identification "CHANNEL" starts in column 1, 
and the actual channel numbers, separated by commas, start in card column 
11, and must be terminated by column 72. 

2. "PRINT" CARD (i„e., "PRINT 1,3,3") 

The "PRINT" card specifies how to print the cluster map. The identification 
"PRINT" starts in column 1, the actual print parameters start in card column 11, 
separated by commss, and ending by column 72. 

3. "HEDl" card 


4. "HED2" card 

These two cards may be used to specify any arbitrary heading for the printer 
output, including the cluster map. Any alphanumeric characters put into 
card columns 11-72 of these two cards will be output as a page heading. 


5. "NPOS" canf 


This card, used previously in the unmodified CLASY program, specified the 
number of positions to skip to read a pixel point (i.e. to "scramble" the 
data). This card is not used in the modified CLASY program, 

6. "NPTS" card 

This card, used previously in the unmodified CLASY program, specified the 
number of pixels to retrieve from the data set at each point. This card 
is not used in the modified CLASY program. (NPTS was used in the original 
programs data scrambling technique). 

7. "DATE" card 

This card is used to specify the date or any eight characters. Will be 
printed at the upper right hand corner of each page of printer output. 

8. "COMMENT" card 

The "COMMENT" card is equivalent in use and format with the "HEDl" and 
"HED2" cards, described above. 

9. "*END*" card 

This card specifies the end of all supervisor (control) card (described 
above) input to CLASY. This card is a mandatory input to CLASY, to 
initiate the clustering process. 

3. 2. 2. 7 Flowchart 



See Appendix A. 


3. 2.2.8 Listings 


See Appendix B for the modified program listing. 

3.2.3 SOFTWARE COMPONENT NO. 3 (CLASYl) 

3. 2. 3.1 Linkage 

CLASYl is the primary entry point, and Is called from CLASY. CLASY2 (an 
alternate entry point In CLASYl) Is called from STATIS. CLASYl calls RREAD, 
CLDUMP, CMERR, and ZOR (Function ZOR Is the random number generator used 
In the data scrambling technique). 

3. 2. 3. 2 Interface 

Interface is accomplished through calling arguments and the following 
common blocks: INFORM, CLUSTR, CLUS, MI SC, and STPAR. 

3. 2. 3. 3 Inputs 


None. 


3. 2. 3. 4 Output 

CLASYl outputs an error message if RREAD attempts to re^iid data from an 
illegal drum address. The following message is printed "***ERROR***RREAD 
TRY TO READ THE ILLEGAL ADDRES', 

3. 2. 3. 5 Storage Requirement 

Storage used: Code = 537g Data = 66g 



3. 2. 3. 6 Description 


CLASYl performs the input image data-handling function for the CLASY clus- 
tering system and makes the image data available to the iterative statistical 
subprogram, STATIS. The original image data from the area on the input tape 
(file) which has been designated by the input field-definition card(s), is 
prestored on drum as one continuous block of data. The Univac random file 
access routines — RINIT, RWRITE — are utilized to place the input data 
on the drum as the data is read from the input tape (file) by the TAPERD 
subprogram. 

CLASYl precomputes the base addresses for three data-buffering arrays in 
core-storage. One data-buffering array, 0, is used to retrieve a block 
of original image data from the prestored drum. 

Another array. A, is used to contain a set of integers — 1,2,3,..., N 
where N=the number of data vectors in the original image array, 0. The 
third array, Dg, will contain a set of data vectors obtained from the 
original image set, 0, but stored in D^, such that each data vector's 
original spatial location is randomly rearranged. The data scrambling 
technique utilized in placing the rearranged data vectors into D^, was ob- 
tained from Rice University via RT&E and is the required modification im- 
plemented by this modification to the CLASY clustering system. The im- 
plementation of the randomized data scrambling technique is performed as 
follows: 

(1) Given an array of original image data vectors, 0, an array, 

A, of integers — 1,2,3, , N with N = the number of data 

vectors in 0, and an array, Dg, for storage of scrambled data 
vectors obtained from 0. 


(2) Scramble the elements of A 


a. Obtain a random number, Zj, from the uniform random 
number generator; 

Zi«Z0R(0), 0.<Z,-<1.0 

b. Multiply the random number Zi, by N, the largest 
integer in A; 

IX^-=N X Z. + 1 

c. Using IX as an index, scramble the integers in A 
as follows: 

TEMP = A^ 

Ai = 

Ajx^ = TEMP 

d. Execute the above procedure (a-c) N times, with 

i = N, N-1, N-2, ..., 1 

e. Create the new (scrambled) set of data vectors, D^, 
as follows; 

for i = 1 , 2 , 3 , . . . , N 

Dc. “ 0(j), where j = Ai 
i 

The scrambled data set, Dg, is made available in large blocks for retrieval 
on demand from the statistical clustering subprogram, STATIC, via the al- 
ternate entry in CLASYl, which is CLASY2. The scrambling of data for STATIS 
is necessitated by the clustering algorithm's sensitivity to correlated data. 

The buffering of blocks of data from the drum to the original data array, 0, 
is done in parallel with the use of scrambled data from by STATIS. The 
Univac random file access routine, RREAD, is used to transfer the data from 
drum to original image array, 0. 


The data buffering technique as described accomplished the second objective 
of the modifications to the CLASY clustering system— namely, the improvement 
of program execution time. 

3.H.3.7 Flowchart 

See Appendix A„ 

3. 2. 3. 8 Listings 

See Appendix B for program listing. 

3,2.4 SOFTWARE COMPONENT NO. 4 (MULTI) 

3. 2.4.1 Linkage 

MULTI is called from CLASY. MULTI calls DATFIX, ALFREE, CLINIT. STATIS 
and CLDUMP. 

3. 2. 4. 2 Interface 

Interface is accomplished through calling arguments and the following 
common blocks: CLUS, MISC, STPAR, INFORM and CLUSTR. 

3. 2. 4. 3 Inputs 
None. 

3. 2. 4. 4 Output 

Upon leaving this subroutine, a message is printed. 


"EXIT FROM MULTI" 


3-2. 4. 5 Storage Requirements 


Storage used: Code = IIIq Data ~ 15j^ 

3. 2. 4. 6 Description 

MULTI calls the routines to initialize the clustering algorithm. MULTI 
was modified to make use of unused areas in core. ARRAY (EXARAY) was pre- 
viously dimensioned but never used. A method was devised to make use of 
this unused core such that the PV array and DATAB array (array utilized in 
STATIS and the array containing scrambled data vectors) be made to utilize 
the same area of core (i.e., the PV and DATAB arrays were made equivalent 
storage areas, in ARRAY (EXARAY)). 

3. 2. 4. 7 Flowchart 


See Appendix A. 

3‘2-4.8 Listings 

See Appendix B for program listing. 

3.2.5 SOFTWARE COMPONENT NO. 5 (STATIS) 

3. 2. 5.1 Linkages 

STATIS is called by MULTI. STATIS calls DISC, CLASY2, CORECT, DOTSQ, 
VPV, VMTV, MPVS, ADJUST, CLDUMP, and EXP. 

3. 2. 5. 2 Interface 


Interface is accomplished through calling arguments and the following 
common blocks: CLUS, MISC, STPAR, CLUSTR, and RAND. 


3. 2. 5. 3 Inputs 


None. 

3. 2. 5. 4 Outputs 

STATIS outputs two warning messages. They are: "***WARNING ON THE__ 

INDEX(KL)=", "**SUSPECTED BAD DATA POINT -STATIS**ID0«__,R00T_, 

VECT0R__" 

3. 2. 5. 5 Storage Requirements 

Storage used: Code = 1205g Data - 260g 

3. 2. 5. 6 Description 

STATIS takes each input data vector and classifies it on a fractional, 
probabilistic basis. It then updates the various statistical parameters 
associated with the classes (clusters) indicated and checks to see if any 
of these classes is potentially two. Those which are will be referred 
to the routine "SPLIT". The one modification made to STATIS was the re- 
dimensioning of the PV array, which contains the data vectors to be clustered. 

3. 2. 5. 7 Flowchart 
See Appendix A. 

3 . 2 . 5 . 8 Listings 


See Appendix B for program listing. 


3.2.6 SOFTWARE COMPONENT NO. 6 (ADJUST) 


3. 2. 6.1 Linkage 

Adjust is called from STATIS. ADJUST calls GET, TR, DOTSQ, SQMTX, MINV, 
UNIF, CLPR, TRIMTX, DENCAL, SPLIT. FREE. CLDUMP, SEPER. SUBLIM, ELIM. 
CORECT, JOIN. APRIOR. SQRT, ALDG, EXP, and XPRI. 


3. 2. 6. 2 Interface 


Interface is accomplished through calling arguments and the following 
common blocks: CLUS, MISC, STPAR, CLUSTR, and JOINPR. 


3. 2. 6. 3 Inputs 


None. 

3. 2. 6. 4 Outputs 

ADJUST prints out three brief messages concerning statistical information 

and three error messages. They are; "ADJUST WEIGHT^ WAS SPFAC 

CHANGE ", "STATISTICS: TRACE_ SKEW__ KURT_ TESTS (SPLIT >0): 

"###HAVE SPLIT_ WEIGHT__ SUBS "W/OVOL ERROR IN ADJUST 

KL.W.NEW W.VOL ", "***EXTRAPOLATION PROBLEM IN ADJUST: ITER, 

INDEX(KL), VOLIN, OVOL, CVOL ", "LOG ERROR IN ADJUST: I, 

IM, KL, K/VRIN= 

3. 2. 6. 5 Storage Requirements 

Storage used: Code = 2305g Data = 332g 


3. 2. 6. 6 Description 


One modificatio made to ADJUST was to eliminate the use of subscripted 
subscripts for those arrays which previously presumed the RFOR (reentrant 
Fortran) compiler. The change was made to enable the routine to compile 
under the Unlyac 1110-EXEC 8 Fortran V compiler. The second change made 
was that Input to the "AL06" routine is forced to be positive by use of the 
absolute value of the Input being sent to "ALOG". The reason for the change 
was that an occasional negative value was‘ being sent to "ALOG", causing an 
error termination. The routine previously referenced a dimensioned variable 
In several calling statements causing an error during program execution. 
These errors were corrected. Other modifications which were suggested by 
the CLASY program originator, Dr. Mike Rassbach, are as follows (where 
"-N" Is Fortran compiler change card format, Indicates that it is a 
change card "N" or "M,N" are the line numbers, in the Fortran symbolic 
listing, to be changed): 

-9 

MQP = MQ+1 
-19 

C Properly, Kurt should be adjusted for the 
C discrete point effect. This has not yet 
C been done, but should not have any 
C major effect, since Kurt Is used 
C only in the crude scan. 

-70, 70 

EXF = WINFC*VACCEL(KADTY) 

-73 

WADJ(KL) = W(KL)*(1.+DWFAC) 

IF(W(KL).LT.WSIM) WADJ(KL) = 2.*W(KL)+WDELSM 

DC0RR=(DW+WADJ[KL])/24 

IF(KADTY.EQ.2) 0C0RR=DW/12. 

DO 118 I = 1, MQS,MQP 
118 ALINK(LB+I-1)=ALINK(LB+I-1) + DCORR 



/c 




-74, 74 

EXF = WINFC * MACCEL(KADTY) 

- 86 , 86 ] 

IF(ITX.GE.10)PRINT 772,ITX,INDEX(KL) ,VOLIN(KL) ,OVOL,CVOL | 

-88 j 

IF(ITX.GE.20)CALL CLPR(KL,NADJ,SUM,SKEW,KURT) j 

-99 

CALL CLPR(KL,NADJ, SUM, SKEW, KURT) 

3. 2. 6. 7 Flowchart 

■ 

I 

See Appendix A. I 

I 

3. 2.6.8 Listings j 

I 

See Appendix B for program listings. | 

i 

3.2.7 SOFTWARE COMPONENT NO. 7 (TR) j 

3. 2. 7.1 Linkage | 

j 

TR is called from ADJUST. | 

1 

I 

3. 2. 7. 2 Interface j 

J 

Interface is accomplished through calling arguments and the following common I 

blocks: CLUS, MISC, and STPAR. ^ 

•j 

I 

3. 2. 7. 3 Inputs 1 


3.2. 7.4 Outputs 


None. 


3. 2. 7. 5 Storage Requirements 
Storage used: Code « 124g Data = 35g 

3. 2. 7. 6 Description 

The modification made to TR was the elimination of subscripted subscripts, 
to enable this routine to be compiled under the Fortran V compiler on the 
Uni vac 1110 (EXEC 8 system). The subscripted subscript notations was a 
feature of the reentrant Fortran compiler RFOR, utilized by the program 
originator in development of CLASY. 

3. 2. 7. 7 Flowchart 


See Appendix A. 

3. 2. 7. 8 Listings 

See Appendix B for program listing. 

3.2.8 SOFTWARE COMPONENT NO. 8 (CLPR) 

3 . 2 . 8 . 1 Linkage 

CLPR is called from CLDUMP, ADJUST, SEPER and JOIN. CLPR calls GET, LOCK, 
SQMTX, MINV, and FREE. 

3. 2. 8. 2 Interface 

Interface is accomplished through calling arguments and the following common 
blocks: CLUS, MISC, and STPAR. 


3 - 


3.2.8. 3 Inputs 


None. 


3. 2. 8. 4 Outputs 

CLPR prints out seven messages concerning statistical information. These 
messages are: "CLUSTER^. INDEX_ PROPORTION^ W*„ SPLIT__ WEIGHT__ WAS_ 

ADJUST__ T0__ PROPORTION: PR0P_ CIN_ CTUT_ OLD PR0P__ CIN_ 0DEN__ 

DIFFER__ VOLUME^, ROUT_ DC0N__", "LOCATION__ LINK__ SUBS__ SUPER__ __ 

SYMB0L__" , "NET PR0B__ DIRECT__ GUMS , "Cl)MS_. , "MEAN , 

"KURT(*W) "OLD COVARIANCE 

3.2.8. 5 Storage Requirements 

Storage used: Code = 740g, Data = 353g 

3. 2. 8.6 Description 

The modifications made to CLPR were (1) the elimination of subscripted 
subscripts (2) a reference to one of CLASY's clustering routines, named 
"LOC", caused ambiguity, because there existed a Univac system routine 
with the same name. Therefore this reference and all other references to 
"LOC" were changed to "LOCK". The originator of CLASY , Dr. Michael 
Rassbach, provided changes to be made in CLPR. One modification suggested 
by Dr. Rassbach was the changing of a format statement so that the printer 
will skip a line before writing the statistical mean. Additionally, the 
following changes were suggested by Dr. Rassbach, and implemented in CLPR: 
-53, 54 

LA=GET(MQS) 

LB-GET(MQS) 

Where "-53, 54" indicates that line No.'s 53, 54 were replaced by the two 
Fortran statements following "-53, 54". 


3. 2. 8. 7 Flowchart 


See Appendix A. 

3. 2. 8. 8 Listings 

See Appendix B for program listing 
3.2.9 SOFTWARE COMPONENT NO. 9 (CLDUMP) 

3.2.9. 1 Linkage 

CLDUMP is called from CLASYl. MULTI. CLASY2 and ADJUST. CLDUMP calls 
ISPLIT and CLPR. 

3. 2. 9. 2 Interface 

Interface is accomplished through calling arguments and the following common 


blocks: 

CLUS, MISC, STPAR, and CLUSTR. 

3. 2. 9. 3 

Inputs 

None. 


3. 2. 9. 4 

Outputs 


A header concerning the dump of the clusters and a debug statement giving 
the value of the number of the split cluster and the value of the print 
parameter. 

3. 2. 9. 5 Storage Requirements 


Data = 33g 



Storage used: Code = 153g, 


3. 2. 9,6 Description 


The modifications made to CLDUMP were the elimination of the logical "IF" 
statement that checked the value of ISPLIT and the print variable, PROUT, 
before the call is made to CLPR. A debug printout statement was added to 
allow printout of these variables for checkout purposes. 

3, 2. 9. 7 Flowchart 


See Appendix A. 

3, 2. 9. 8 Listings 

See Appendix B for program listing. 

3.2.10 SOFTWARE COMPONENT NO. 10 (ELIM) 

3.2.10.1 Linkage 

ELIM is called from ADJUST. ELIM calls SUBLIM and TRFREE. 

3.2.10.2 Interface 

Interface Is accomplished through calling arguments and the following 
common blocks: CLUS, MISC, and STPAR. 

3.2.10.3 Inputs 
None. 



3.2.10.4 Outputs 


ELIM prints out two messages, one statement when a cluster has been eliminated 
and the other is an error message when a structural error has occurred. 

These messages are: ''###ELIMINATE__ LINK, LSUBS, LSUPER 

"♦♦STRUCTURAL ERROR AT ELIM; KEL, KFAITH, KOLD, INIT 

3.2.10.5 Storage Requirements 

Storage used: Code = 176g, Data « 46g 

3.2.10.6 Description 

The modifications made to ELIM were the eliminations of subscripted sub- 
scripts in arrays to allow the routine to be compiled under the Fortran V 
compiler of the Univac 1110-EXEC 8 system. 

3.2.10.7 Flowchart 


See Appendix A. 

3.2.10.8 Listi nqs 

See Appendix B for program listings. 

3.2.11 SOFTWARE COMPONENT NO. 11 (SEPER) 

3.2.11.1 Linkages 

SEPER is called from ADJUST. SEPER calls CLPR, DGNCAL, and FREE. 



3.2.11.2 Interface 


Interface is accomplished through the following common blocks: CLUS, MISC, 
and STPAR. 


3.2.11.3 Inputs 
None. 

3.2.11.4 Outputs 

SEPER prints one statement, which is printed whenever a cluster Is split. 

This statement is ; "###SEPARATE_ SUPER. SUBS SPFAC_'' 

3.2.11.5 Storage Requirements 
Storage used; Code = 242g, Data * 46g 

3.2.11.6 Descriptions 

The modifications made to SEPER were the elimination of subscripted sjbscripts 
in arrays, to enable the routine to be compiled under the Fortran V compiler 
of the Uni vac 1110-EXEC 8 system. 

3.2.11.7 Flowchart 
See Appendix A. 

3.2.11.8 Listings 



See Appendix B on program listing. 


3.2.12 SOFTWARE COMPONENT NO. 12 (LOCK, formerly “LOC) 


3.2.12.1 Linkages 

LOCK is a function subprogram and is called from CLPR, and CLPRM. 

3.2.12.2 Interface 

Interface is accomplished through the function arguments. 

3.2.12.3 Inputs 
None. 


3.2.12.4 Output 


None. 


3.2.12.5 Storage Requirements 
Storage used: Code == 16 ^q 

3.2.12.6 Descriptions 

The modification made to LOCK (formerly "LOC") was the changing of the 
external reference name from 'LOC to 'LOCK'. The name 'LOC is also the 
name of a Fortran V system routine which caused ambiguity between the two 
routines during execution of CLASY. Changes were made to all other routines 
that referenced this routine, to change the reference to be "LOCK" rather 
than "LOC". 

3.2.12.7 Flowchart 


'lU 


See Appendix A. 


3.2.12.8 Listings 


See Appendix B for program listing. 
3.2.13 SOFTWARE COMPONENT NO. 13 (MISH) 

3.2.13.1 Linkages 


The Fortran V procedure, MISH, is included in the following routines: 
ACOM, ADJUST, ALFREE, AMSQ, APRIOR, CBLO, CLASY, CLASYl , CLDUMP, CLINIT, 
CLPR, CLPRM, CLUSMP, CLUST, CORECT, DATFIX, DOTSQ, EIGROT, ELIM, ISPLIT^ 
JOIN, MINV, MLT, MPVS, MTVEC, MULTI, MVEC, SEPER, SPLIT, SQMTX, STATIS, 
STOFLO, SUBLIM, TR, TRFREE, TRIMTX, VMTV, and VPV. 


3.2.13.2 Interface 


The Fortran V procedure, MISH, is used as the interface for the following 
common blocks: CLUS, MISC, and STPAR. 


3.2.13.3 Inputs 


None. 

3.2.13.4 Output 
None. 


3.2.13.5 Storage Requirements 


Storage used: Code ~ 50 


3.2.13.6 Descriptions 


The modification made to the proc, MISH, was the changing of the symbol 
array *NSYMB' dimension from 11 to "12" (i.e., NSYMB(12)). This was done 
because another symbol was added to the large symbol array, *SYM'. 

3.2.13.7 Flowchart 


See Appendix A. 

3.2.13.8 Listings 

See Appendix B for program listing. 

3.2.14 SOFTWARE COMPONENT NO. 14 (CMBKIO) 

3.2.14.1 Linkages 

The Fortran V procedure, CMBKIO, is included in the following routines: 
ADJUST, CLASY, CLASYl , CLDUMP, CLINIT, CLUSMP, JOIN, MULTI, READTP, SETUP9, 
and STATIS. 

3.2.14.2 Interface 


The Fortran V procedure, CMBKIO, is used as the interface for the common 
block , CLUSTR. 

3.2.14.3 Inputs 

None. 


3.2.14.4 



None. 


3.2.14.5 Storage Requirement 


Storage used: Code “ 10 

3.2.14.6 Description 

The modification made to the Fortran V procedure CMBKIO was the changing of 
the parameter, MAXPOP, value from 60 to 61. This was done to increase the 
dimension of the symbol array to account for the addition of another symbol 
to the array. The symbol added to the array was the character 'O' (zero). 

3.2.14.7 Flowchart 
See Appendix A. 

3.2.14.8 Listings 

See Appendix B for program listing. 

3.2.15 SOFTWARE COMPONENT NO. 15 (CLUST) 

3.2.15.1 Linkage 

CLUST is called from CLUSMP. CLUST calls ISPLIT, CORECT, DOTSQ, and EXP. 

3.2.15.2 Interface 


Interface is accomplished through calling arguments and the following common 
blocks: CLUS, MISC, STPAR, BIGCOM. 

3.2.15.3 Inputs 

None. 

t 
( 

t 

i 

( 

} 



3.2.15.4 Outputs 


None. 

3.2.15.5 Storage Requirements 

Storage used: Code » 225„, Data = 56 q 

8 o 

3.2.15.6 Descriptions 

The originator of the CLASY clustering system, Dr. Michael Rassbach, pro- 
vided the following changes to incorporate into the routine 'CLUST'. The 
changes suggested by Dr. Rassbach, and implemented in CLUST are as follows: 
-15 

C USES PCOND FOR PPASS, DIST FOR DISS (KL) 

-18, 18 

PCOND (KROT) = 1. 

-25, 25 

130 PCOND (KL) »= PROP (KL)/PR IRCM (KFAITH)*PCOND(KFAITH) 

-34,34 

DIST = DOTSQ (REL,VRIN(KL+1)*W(KD) 

-35,35 

IF(ABS(DIST+DC0N(KL)).LE.160.) GO TO 531 
-38,38 

P = XP(DIST+DCON(KL))/VOLRT(KL))*PCOND(KL) 

-46 

149 PCOND(KL) = 0 
-47,47 

KL = KFATH 



3.2.15.7 Flowchart 


See Appendix A. 

3.2.15.8 Listing 

See Appendix B for program listing. 


4. OPERATION 


4.1 TEST PROCEDURE 


The test procedure for verification of the modifications to the CLASY clus 
tering system is described below. 

For the identical image data set and CLASY input run parameters (area of 
image clustered, and channels selected): 

(1) Provide a baseline run on the original (pre-modification) CLASY 
clustering program, 

(2) Provide a second run on the CLASY clustering program with all 
modifications documented herein, 

(2a) Compare program execulion times for the two runs, for 
an indication of improvement. 

(3) Provide a run on the modified CLASY program which will list the 
array of input original image data, and the array containing the 
input data to be passed to STATIS after application of the 
randomized data scrambling technique, 

(3a) Compare the listing of the two arrays for any evidence of 
data value changes or extraneous data values introduced. 

(4) Compare the cluster map output from the modified CLASY program to 
thfe cluster map output by the baseline (original) program, with re- 
ference to the known statistics and class definitions of the simu- 
lated image data used to make the runs. 

4.2 CLASY PROGRAM OPERATION 

The CLASY clustering program, as modified, is operational on the Univac 
1110 (or 1108) under the EXEC8 operating system. The program utilizes 
the Univac Fortran V compiler (the original program utilizes the Univac 


reentrant Fortran compiler, RFOR), the Univac assembler (for three assembly 
language routines—FREE, GET, and LOCK), and the Univac system random file 
access routines, RINIT, RREAD, and RWRITE. 

4.2.1 CLASY Program Input 

The required input to the CLASY program consists of one tape (or file) con- 
taining the multichannel image data, and the special-format card input. 

The image data tape (file) is presumed to be in either of two specific 
formats--etther "LARSYS II" format or "UNIVERSAL" format. The tape (file) 
reading program in CLASY, TAPERD, accepts either of these formats and self- 
determines the correct method of reading the data. 

The card input to CLASY is of a specialized "free-field" format. Two types 
of card input are accepted--"control" cards (not to be confused with Univac 
system control cards) and "field definition" cards. The control cards are 
used to implement the execution options available in the program. The key- 
words for each "control" card are shown below. The field definition cards 
are used to delineate the specific areas of the input image data which are 
to be clustered by CLASY. The format of the field definition cards is 
shown below. 

4. 2. 1.1 Control Card Input 

The format and use of the "control" cards is as follows: 

(1) The keyword begins in card column 1, and must be spelled exactly 
as shown. 

(2) The parameter(s) on the control card begin in or after card column 
11 and end in or before card column 72. 

(3) The comma is the recognized separator for integer parameters, 
and blanks are ignored. 

(4) If a default is shown, the control card is optional input to the 
program (default value is provided if the card is not input). 



(5) If the default is specified as "none", the card is mandatory input to 
the program. 

(6) The order of the input control cards is immaterial, except for the 
"*END*" card and the "$END*" c»>rd. The "*END*" card must be the 
last card of the set of control cards, and is to be followed im- 
mediately by one or more field definition cards. The "$END*" card 
must follow the last field definition card. The "$END*" card is the 
last card of the input. 


(7) The "control" cards recognized by CLASY are as follows: 


\ # ... — 

Keyword 

Paramater(s) 

Function 

(CC 1-8) 

(CC 11-72) 


HEDl 

Any alphanumeric 
characters (default: 
standard page heading, 
line 1) 

The characters supplied in 
CC 11-72 will replace the 
first line of output page 
heading 

HED2 

Any alphanumeric char- 
acters (default: stan- 
dard page heading, line 2) 

The characters supplied in 
CC 11-72 will replace the 
second line of the standard 
output page heading. 

DATE 

Any eight alphanumeric 
characters (default: 
the current calendar 
date, in upper right-hand 
corner of each output page) 

The (up to) eight characters 
beginning in CC 11 will be 
placed in the upper right-hand 
corner of each page of printer 
output. 

CHANNELS 

Nl, N2, N3, N4 with N an 
integer (default: none) 

The channel numbers to be used 
in selecting the data vectors 
to be clustered. The channel 
numbers must be a subset of the 
actual channels available in 
the input image data. 

PRINT 

1,3,3 

(default: limited in- 
terim printout during 
program execution) 

The parameters "1", "3", 

"3" are used as printout con- 
trol to generate diagnostic 
data during the iterative 
clustering process. 


Keyword (Cont*d) Parameter(s) (Cont'd) Function (Cont*d) 


*END* None 

(default: none) 


SEND* None 

(default: none) 


Mandatory card-- signifies the 
end of the control card set, 
and the beginning of the field 
definition card input. 

Mandatory card- -must follow 
the last field definition 
card— initiates CLASY program 
execution. 


4. 2. 1.2 Field Definition Card 

The "field definition" card(s) delineate the area on the image data tape 
(file) which is to be clustered by CLASY, in terms of pixel coordinates 
(sample, line) for each vertex of the "field", or area, up to a maximum of 
10 vertices for a given field. An alphanumeric field identification may be 
supplied in card columns 1-6, but is not required. Beginning in card 
column 11, through card column 72, coordinate pairs are given (always in 
sample, line order) enclosed in parentheses, separated by a comma The 
first pair given for a field must be the incrementation desired in the 
lines and pixels to be read from the input image tape (file). I.E., "(2,3)" 
would indicate every second pixel on each line, and every third line to be 
read. The second and succeeding coordinate pairs are the (sample, line) 
coordinates of the vertices of the field. A continuation of coordinate 
pairs on the next card is indicated by an asterisk Up to 10 coordinate 

pairs (vertices) will be accepted for one field. An example of the field 
definition card(s) is as follows: 

CC 1 CC 11 

CARD 1 FLDNAM (2,3) , (2,1) , (196,3) , (100,50) * 

CARD 2 (196,100) , (1,100) , (20,30) 

The result of the above card will be that a rectangular area will be read 
from the image data which bounds the given irregularly shaped field defined 
above. The coordinates of the rectangular area will be (1,1) , (196,1) , 
(196,100) , (1,100). 


The coordinates for the rectangular area are self-detennined by the tape 
(file) reading program. Internally, the actual field coordinates which 
were input are used to extract for clustering only the specific pixels that 
are within the actual field defined. 

The input field vertices must be defined on the card(s) in clockwise order . 

4.2.2 CLASY PROGRAM OUTPUT 

The output by CLASY is all 1 ine-printer output. Interim printout of sta- 
tistical parameters and diagnostic data is provided during the iterative 
cl us ter- forming process. 

The final output is a "map", with a symbolic representation of the area 
clustered, with each pixel of the area classified using the statistics 
(mean and covariance) from final cluster set determined by CLASY. The 
symbols on the "map" represent the cluster {- class) which is the most 
likely parent distribution for the given pixel. The "map" is output by 
subprogram CLUSMP. 

Sample output is shown in Appendix C. 

4.2.3 CLASY PROGRAM RESTRICTIONS 

The known restrictions inherent in the program are (1) the program is re- 
stricted to a maximum of four channels, due to program array dimensioning 
which presumes this maximum, (2) the program will not successfully execute 
with only one channel, (3) a data vector containing a zero value in the 
channel of interest will cause an error termination of the program's execu- 
tion, (4) the size of the original image data set read from the input 
tape (or file) and placed on drum must be containable in 1,310,717 locations 
of drum storage available to the random access routines (RINIT, RREAD, 
and RWRITE). 


TEST VERIFICATION 


Po** CLASY Program Modification 


This verification Is being conducted to Insure that the delivered 
program products satisfy the requirements as originally stated by 
the requesting organization. 





Cognizant System Manager 
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